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Abstract 

We report recent results on the properties of the X(3872) produced via the B + — ► K + X(3872) 
decay process. We observe decays X — ► ir + ir~ it J / 'ip where the 3tt invariant masses cluster near 
the upper kinematic boundary suggesting that they originate from sub-threshold decays to virtual 
a; (782) mesons. This is consistent with expectations for a DD* bound state interpretation for the 
X(3872). In addition, we constrain the possible charmonium-state assignments for this particles. 
Results are obtained from a 253 fb -1 data sample that contains 274 million BB pairs that was 
collected near the T(4S) resonance with the Belle detector at the KEKB asymmetric energy e + e~ 
collider. 

PACS numbers: 



3 



INTRODUCTION 



A first step in understanding the X(3872) particle that was recently discovered by Belle T] 
is to determine its J PC quantum numbers. Here, we survey possible assignments and proper- 
ties of the most likely candidates and contrast these with recent experimental measurements. 

Although the AT (3872) is above the DD mass threshold, its width is narrow, V < 
2.3 MeV pj, and decays to DD are not seen 0. This suggests that DD decays are forbid- 
den. We restrict our considerations to ++ and 1 -k + -k~ systems [3j and scenarios where 
the relative orbital angular momentum of the 7r + 7r~ and J/ip is L < 3. We concentrate on 
possible charmonium assignments, and only those where decays to DD are forbidden or ex- 
pected to be strongly suppressed. For the case of a ++ dipion, there are three charmonium 
possibilities: the /i c (2 1 P 1 ) and two triplet D-wave states, the ip2^D 2 ) and ip^D^). For the 
1 dipion case, there are also three possibilities: the r\ c , the x c i an d the 77 c2 ( 1 Z? 2 ) - For these 
assignments, the 7i + 7i~J/ip decay would violate isospin and should be strongly suppressed. 

The Belle experiment observes B mesons produced by the KEKB asymmetric energy 
e + e~ collider [3]. KEKB operates at the T(45) resonance (y/s = 10.58 GeV) with a peak 
luminosity of 1.39 x 10 34 cm _2 s _1 . At the T(AS), BB meson pairs are produced with no 
accompanying particles. As a result, the B mesons have a total center-of-mass system (cms) 
energy that is equal to Ebeam, the cms beam energy. We identify B mesons using the beam- 
constrained mass M bc = J E^ eam — p B and the energy difference AE = E beam — E B , where 
Pb is the vector sum of the cms momenta of the B meson decay products and E B is their 
cms energy sum. The experimental resolution of M bc is approximately 3 MeV; that for AE 
is typically 11 MeV for all-charged-particle final states. For final states with 7's or 7r°'s, 
the AE resolution becomes broader and somewhat skewed to negative values due to energy 
leakage out of the back of electromagnetic calorimeter. The Belle detector is described in 
ref. 



SEARCH FOR X(3872) -> 7 x c2 (Xcl) 

The Wigner-Eckart theorem says that the widths for ip2 7i + 7i^J/ijj and ip3 — *> tt + tt~J/ 
should both equal r(?/>(3770) — > 7i + 7i^J/tp). This has been recently measured by BESII 
and CLEO-c 7] to be 80 ± 32 ± 21 keV and < 55 keV (90% CL), respectively. The results 
are in some contradiction with each other. For the following discussion we conservatively 
assume an upper limit derived from the larger BES number of r(?/>(3770) — > 7i + 7i^J/ip) < 
130 keV. 

Calculations of the 7Xci width for an M = 3872 MeV ip 2 range from 207 keV to 
360 keV [9]. The 90% CL upper limit of 

that was reported in ref. [1] contradicts these expectations for the f/V 

Barnes and Godfrey |9j observe that although ^3 — > DD is allowed for a 3872 MeV ips, 
this mode is suppressed by an L = 3 centrifugal barrier and the total ^3 width may be less 
than the T < 2.3 MeV experimental upper limit. These authors, and also Eichten, Lane and 
Quigg 0, propose the ^3 as a charmonium candidate for the V(3872). 

For an M = 3872 MeV ^3, the calculated 7% C 2 widths range from 299 keV[8] to 
370 keV 0. Thus, the partial width for ip 3 — * 7x c2 is expected to be more than twice 
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that for -?/>3 — > 7T + 7T~J/ip. We performed a search for X — > 7Xc2 that followed closely the 
procedure used for the 7Xei limit reported in ref. [1]. We used a 140 fb _1 data sample, which 
contains 152 million BB pairs. We require one of the 'yJ/ip combinations to satisfy 444 MeV 
< (M je + t - - M e+e -) < 469 MeV. The M bc and AE signal regions are \M bc - 5.28| < 0.009 
GeV and -0.04 < AE < 0.03 GeV. 

We use the B — > Kip ; ip — > 7Xc2 decay chain as a normalization reaction. The signal- 
band projections of M bc and M 7Xc2 for the ip region are shown in Figs. [T] (left) and (right), 
respectively, together with curves that show the results of the fit. The fitted signal yield is 
18.3 ± 5.2 events, where, based on known branching fractions, we expect 12 ± 3 events. 

Figure El show the same projections for events in the X(3872) mass region, where there 
is no apparent signal. An unbinned fit produces a signal yield of 2.9 ± 3.0 ±1.5 events, 
where the first error is statistical and the second systematic. The latter is estimated by the 
changes that occur when the input parameters to the fit are varied over their allowed range 
of values. 

The ratio of the X —>■ r )Xc2 and the X — > 7t + tt~ J /if) partial widths and its 90% CL upper 
limit are 



0.42 ± 0.45 ± 0.23 < 1.1(90%CL), 



(2) 



where the second quoted error is the quadratic sum of the systematic uncertainties in ac- 
ceptance, the branching fractions and variations in the 7Xc2 event yield for different fitting 
methods. 
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FIG. 1: Signal-band projections of M bc (left) and M 7Xc2 (right) for events in the ip region with 
the results of the unbinned fit superimposed. 



SEARCH FOR X -► 7 J/ip 

The x'ci i s expected to be near 3968 MeV, well above the DD* threshold, and its width is 
expected to be hundreds of MeV@|. If potential models are wrong and the x' c i * s below the 
DD threshold at 3872 MeV, it could conceivably be narrow and 7r + 7r~ J/ 'if) decays might 
be significant, even though these would violate isospin. In this case, the jif) and •yJ/ip 
transitions would be important and almost certainly have larger partial widths than that 
for the 7r + 7r~ J I if) mode. We searched for the X — > jJ/if) decay mode. 

We select B + — > K +r yJ/if) event candidates using the criteria given in ref. [1]. The 
B + — > K +r yJ/ip channel is dominated by B + — > K + x c i, Xci — > iJ/^ decays and we use this 
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FIG. 2: Signal-band projections of M\, c (left) and M 1Xc2 (right) for events in the X(3872) region 
with the results of the unbinned fit superimposed. 




FIG. 3: Signal-band projections of Mf, c (left), M^ji^ (center) and AE (right) for events in the 
Xci region with the results of the unbinned fit superimposed. 

as a calibration reaction. We define a Xci window for jJ/i/j masses within 20 MeV of the 
nominal Xci mass. Figure HH shows the signal-band projections for M bc (left), M 7 jm (center) 
and AE (right) for events in the Xci region with the results of a three-dimensional unbinned 
fit superimposed. The fitted number of events is 470 ± 24. 

We define an X -> jj/ip signal region to be |M (7J/VO - 3872 MeV| < 26 MeV. FigureH 
shows the same projections for events in the X(3872) signal region. Here there is no strong 
evidence for a signal: the fit gives a 2.2<r signal yield of 7.7 ± 3.6 events. The resulting limit 
is 

hA ~ : ' //!i - 0.22 ±0.12 ±0.06 < 0.40(90%CL), (3) 



T(X -> 7T + 7T- J/lp) 



where the second quoted error is systematic and includes uncertainties in acceptance, the 
branching fractions and variations in the jJ/ip event yield for different fitting methods. 




FIG. 4: Signal-band projections of M bc (left), M^ji^ (center) and AE (right) for events in the 
X(3872) signal region with the results of the unbinned fit superimposed. 
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HELICITY ANGLE DISTRIBUTION FOR 1+ ~h' c 



The h' c hypothesis makes reasonably specific predictions for the X — > 7i + 7i^J/ip decay 
angular distributions (loj . We define 9j/^ as the angle between the J/ip and the negative 
of the K + momentum vectors in the X(3872) rest frame in the decay B — > XK; X — > 
it + it~ J/ip. The | cos9j/^\ distribution for X(3872) events with m^^- > 0.65 GeV is shown 
as data points in Fig. El The smooth dotted curve is a polynomial represents sideband- 
determined backgrounds. 




0.00 0.25 0.50 0.75 1.00 

|cose J)¥ i 



FIG. 5: The measured | cos8j/^\ distribution. The superimposed histogram is the normalized MC 
distribution for the 1+- hypothesis. Here x 2 /dof = 75/9. 

Figure shows a comparison of the measured \cos9j/^\ distribution with a MC sample 
generated with a J PC = l + ~ hypothesis. Here the expected | cos 9j/^\ distribution has a 
sin 2 #j/^ dependence that goes to zero at cos9j/^ = 1, where the data tend to peak. This 
makes the overall x 2 quite poor, x 2 /dof is 75/9, and enables us to rule out the l + ~(h' c ) 
hypothesis for the X(3872) with high confidence. 



SEARCH FOR V(3872) ir°ir°J/ip 

The ratio TZq = T(X — > i[ iT Q J/ip)/Y(X — > n + ii^ J/ip) measures the isospin of the dipion 
system If isospin is conserved, 1 = corresponds to 7Z = 1/2; for / = 1, 1Z = 0. For 
the ijj', where the dipion system is known to have / = 0, this ratio is 0.60 ± 0.05 We 



searched for X — > 7r 7r J/?/> decays using a 253 fb~ data sample that contains 274 million 
BB pairs. We use B — > Kip'; %p' — > 7i°TT°J/ifj as a calibration reaction. 

We select B — > Kir 7r J/ip events using the J/ip and charged kaon criteria given in ref. [1]. 
For neutral Bs we use the standard Belle K$ — > tt + tt~ criteria and require the 7t + tt~ invariant 
mass to be within ±15 MeV (~ 3a) of M^o. We identify 7r°s as 77 pairs that fit the ir° —> 77 
hypothesis with x 2 < 6. We further require the energy asymmetry \(E J1 —E 72 )/(E rL +E ry2 )\ < 
0.9 and the 7r° cms momentum to be greater than 150 MeV. For cases where there are more 
than two 7r° candidates in the \ AE\ < 0.2 GeV and M bc > 5.2 GeV selection region, we take 
the combination with the minimum value of (M 77 - M w o) 2 /a 2 + \AE\ 2 /a 2 AE . 
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FIG. 6: (a) The M& c distribution for candidate B — ► i^^'; ^' — ► ir°TT°J/ip decays. The curve is 
the result of the fit described in the text, (b) The Mb c distribution for B — > KX; X — > 7r°7r°J/V' 
decays. 



Figure Ufa) shows the M& c distribution for events in the —0.06 GeV < AE < 0.03 GeV 
signal region with M^o^oj/^ within ±15 MeV of M^>. Here we have also required M^o^o > 
390 MeV. We fit this distribution with a Gaussian function to represent the signal and a 
smoothed threshold function to represent the background. The signal yield from the fit is 
55 ± 10 events. 

Figure Efb) shows the corresponding M bc distribution for events with M^oj/^ within 
±15 MeV of 3872 MeV. We have applied an additional requirement M n o^o > 570 MeV. This 
dipion mass restriction has ~100% acceptance for X — > 7i + 7i^J/ifj decays and we assume a 
similar efficiency for 7r 7r J /ip. Here there is no evident signal; a fit gives a signal yield of 
0.2 ± 2.6 events. 

We compare these signal yields to corresponding results for ip' and V(3872) decays to 
7r + 7r~ J/ip and determine the 90% CL upper limit: 

By quoting our result this way, many systematic effects, including the influence of the dipion 
mass requirements, cancel out. Unfortunately, with the present data sample, the limit is 
not stringent enough to distinguish between the 1 = and 1 = 1 hypotheses. 



OBSERVATION OF V(3872) -> ujj/ip 

In the context of an analysis of the possibility that the X(3872) might be a weakly 
bound J PC = 1 ++ D°D°* molecular state, Swanson [lj| developed a specific model for the 
X(3872) as a D°D°* hadronic resonance with an important admixture of uJ/ip. He finds 
that although the X is 7.5 MeV below the Mj/^ ± M u mass threshold, decays that proceed 
via virtual cu mesons are important; the signature being B — > KX; X — > 7r + 7r~7r° <//■?/>, where 
all the events have 7r + 7r~7r° invariant mass values near the kinematic upper limit. He predicts 
that these decays should occur at about half the rate for X — > 7r + 7r~ J/ip. We searched for 
B — > KX\ X — > 7r + 7r~7T° J/i/j decays in a 274 million BB event sample. 
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We select events with the same kaon, charged pion and J/ip requirements as used in 
ref. [1]. An additional 7r° — > 77 candidate that satisfies the same criteria listed above for the 
X — > 7t TT°J/ip search is required. The tt° cms momentum is required to be above 180 MeV. 
We select events where M(n + n^n°J/ip) is within ±12 MeV (2a) of 3872 MeV. The M bc and 
AE signal regions are \M bc - M B \ < 7.5 MeV and \AE\ < 30 MeV. 




FIG. 7: (a) The M bc and (b) AE distributions for candidate B -> K-k + tt~tt q J decays. The 
curves are the result of the fit described in the text, (c) The M(7r + 7r~7r°) distribution for events 
in the M bc -AE signal region. 

Figure Ufa) shows the M bc projection for events in the AE signal region; Fig. Efb) shows 
the corresponding AE projection. There is some indication of a B meson signal: a simulta- 
neous fit to the two distributions gives a signal yield of 15.4 ± 6.3 events, with a S/B — 0.8. 
Figure UJc) shows the 7T + 7r~7r° invariant mass for events in the M bc -AE signal region, where 
a peak is evident at high masses. 
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FIG. 8: (a) The M bc and (b) AE distributions for candidate B — > Kir + it ir°J/ip decays with 
M(7T + 7r~7r°) > 0.75 GeV. The curves are the result of the fit described in the text. 

Figures Ufa) and (b) show the M bc and AE projections for events with M(7r + 7r~7r°) > 
0.75 GeV (indicated by the arrow in Fig. [7(c)). Here a B meson signal is evident on a small 
background. The signal yield from a simultaneous fit is 10.0 ±3.6 events and S/B = 5. The 

statistical significance of the signal, determined from ^J— 2 ln(£ / Cm)) where £ max and £ 
are the likelihood values for the best-fit and for zero-signal-yield, respectively, is 5.8cr. This 
is the first obervation of an V(3872) decay mode other than 7i + 7i^J/ip. 

We attribute all of the signal events with 7r + 7r~7r° invariant mass greater than 0.75 GeV 
to B — > ujJ/ip and compute the ratio of uuJ/ip and tt + tt^J/iIj partial widths by comparing 



9 



this to the number of X — > 7r + 7r J /if) in the sample data sample, corrected by the relative 
detection efficiencies: 

TiX^uJ/ib) . 

0.8±0.3(stat)±0.1(syst), (5) 



where the systematic error reflects the uncertainty in the relative acceptance. 

The properties of the X — > 7r + 7r~7r° J/^ decays are in good agreement with expectations 
of ref. |l3j. The 7r + 7r _ 7r° invariant masses cluster near the upper kinematic limit, and its 
measured strength is consistent with being "roughly 1/2" that for the tt + tt~ 'J /if} mode. 



SUMMARY 

The observation of X(3872) — > n + n^n J/ip decays with properties consistent with expec- 
tations for X — > uJ /if) provides strong support for the DD* molecular state interpretation 
for the X(3872). Both the observed 37r mass distribution and the decay srength relative to 
n + n^J/ip agree with predictions by Swanson 13]. 

Moreover, none of the six possible charmonium candidate states comfortably fit the mea- 
sured properties. The 90% CL branching fraction upper limit for B(X (3872) — > 7X c s) 
decay is 1.1 times that for 7r + 7r~ J/-?/>. This conflicts with theoretical expectations for the 
case where the X(3872) is the 3 ips- 

The possibility that the X(3872) is the 1 ++ X c i charmonium state is made improbable by 
the limit B(X — > jJ/if>) < 0AB(X — > 7r + 7r~ J/ifA- The former would be an allowed El tran- 
sition with an expected width of T lJ /^ ~ 10 keV|9j. The latter would be an isospin- violating 
transition; other isospin violating transitions in the charmonium system have widths that 
are less than 1 keV. 

An analysis of the 6j/^ helicity angle distribution eliminates the l + ~(h c ) hypothesis with 
a high degree of confidence. 

The 0"+(^) mass differs from that of the ip(3S) by hyperfine splitting and can be reliably 
expected to be about 50 MeV (or less) below that of the if>(3S), which is at 4030 MeV. 
Moreover, even if it were as low as 3872 MeV, the width is expected to be some 10 's of 
MeV, similar to that of the r\ c and wider than the 2.3 MeV upper limit for the X(3872). For 
2~ + (?7c2), the T] C 2 — > n + 7i~r] c and jh c decays are allowed and expected to have widths in the 
range of 100's of keV[9], and much larger than that for the isospin- violating 7r + 7T~J/ip mode. 
If the X(3872) were the i] C 2, the total exclusive branching fraction for the B + — > K + r] c2 decay, 
which is non-factorizable and suppressed by an L = 2 barrier, would be anomalously large. 
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